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We report preliminary results of a Wilson Chiral Perturbation Theory (WChPT) analysis of 
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was generated with two definitions for the critical quark mass and shows a strong non-linear quark 
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estimates for combinations of low-energy parameters, even though the errors are sizable. 
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1. Introduction 

Twisted mass Lattice QCD (tmLQCD) [1,2] is a promising formulation to approach the chiral 
limit of QCD with Wilson fermions. It solves the problem of exceptional configurations and makes 
numerical simulations with small quark masses feasible. Recent quenched studies [3-6] reached 
values for m n /m p as small as 0.3 without running into problems due to exceptional configurations. 
Besides this numerical advantage tmLQCD has the property of automatic 0(a) improvement at 
maximal twist [7]. More recent results in tmLQCD have been reviewed at this conference [8]. 

Some issues, however, remain to be fully understood. In order to be at maximal twist the 
bare untwisted mass mo needs to be tuned properly. It has been argued in Ref. [9] that automatic 
0(a) improvement is lost if one tunes to the critical quark mass where the pion mass vanishes, 
unless the twisted quark mass satisfies the bound > a 2 AQ CD . On the other hand, automatic 
0(a) improvement is expected to hold if mo is tuned such that the PCAC quark mass vanishes. 
However, the results in Ref. [9] are controversial. It is claimed [6, 10] that both definitions for the 
critical quark mass guarantee automatic 0(a) improvement, but the remaining 0(a 2 ) effects differ 
significantly. In particular, the pion mass definition is afflicted with cut-off artifacts of 0(a 2 /m 2 ) 
which are enhanced for small pion masses [10]. These enhanced lattices artifacts are shown to be 
reduced for the PCAC quark mass definition. 

Closely related is the so-called 'bending phenomenon', observed in quenched simulations 
[3-6]. Various mesonic quantities like the pion mass, the pion decay constant and the vector meson 
mass show an unexpected non-linear quark mass dependence for small quark masses if the pion 
mass definition for the critical quark mass is used. This curvature is absent when the PCAC quark 
mass definition is employed. 

Here we report preliminary results of a Wilson Chiral Perturbation Theory (WChPT) analysis 
of the lattice data in Refs. [3-6]. Our motivation for this analysis is two-fold. Firstly, we would like 
to check whether WChPT is able to describe the lattice data and the peculiar bending, as it should 
if it is the correct low-energy theory of tmLQCD. Previous analyses with WChPT [1 1-13], applied 
to data obtained with untwisted Wilson fermions, are inconclusive and do not provide sufficient 
evidence for the validity of WChPT. Secondly, provided WChPT can be applied, a fit to the data 
should give us some combinations of low-energy coefficients of WChPT. Particularly interesting 
is the coefficient C2 which plays a crucial role for the phase structure of the lattice theory [14-16]. 
This coefficient can be determined by measuring the mass difference between the charged and 
neutral pion [17], however, this calculation is difficult in practice since it involves the computation 
of disconnected quark diagrams. A recent computation of the pion mass difference [18] shows 
large statistical errors and alternative methods to get a handle on C2 are clearly welcome. 

2. WChPT for twisted mass LQCD 

We extended the WChPT analysis in Ref. [9] by including the terms of 0(/xa) and 0(a 3 ) in 
the chiral Lagrangian. Keeping these terms is consistent in a power counting scheme where the 
renormalized masses m and fl and the lattice spacing a satisfy (proper powers of Aqcd are , as 
usual, understood) 

jj. ps a 2 > \la ps a 3 > \l 2 ps jla 2 ss a 4 . (2.1) 
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Even though the 0(pia,a 3 ) contributions are subleading, we found them to be essential to describe 
some features in the lattices data. 

A twisted mass term breaks the flavor symmetry. The ground state £o of the theory is therefore 
no longer symmetric under flavor rotations but given by 

I = exp(/0T 3 ), (2.2) 

where denotes the vacuum angle. Minimizing the potential energy we derive a gap equation, 

a 2 (?-/3(l-2? 2 )) 2 = (x-t+iPXnf + ytfiX-t 1 ), (2-3) 

which determines t = cos0 as a function of pi, m and a. The masses and the lattice spacing are 
included in the dimensionless parameters 

2Bpi 2Bm + 2Woa(l +cj,a 2 ) 2Bm c%a 3c3« 3 

^ = o ?i X = ~ o ) Xm = ~ ?) P = 7 = o t")(2-4) 

2cia L 2c2a z 2c2a z 2B 2c2<2 

where B, Wq, c ; -, c,- are (combinations of) low-energy parameters in the chiral Lagrangian. Expand- 
ing in terms of pion fields in the usual way one derives the pion mass, the pion decay constant and 
the twist angle odwt stemming from the chiral Ward identities [1]: 



2 2B\x 1 + fit 



1-Sln 



25/1 



(2.5) 



/* = /(! + S[t + S^y/l^fi) Vl^, (2-6) 

cotaiwT = , = , (2.7) 

VI -f 2 

where 8\ , 5( , 5^ are additional combinations of low-energy coefficients and 8 multiplies the quenched 
chiral log contribution [19,20]. For f K we also included the 0(/i 2 ) correction which we dropped in 
the result for the pion mass. 



The results ( |2.5D - ( |2.7| ) are valid for arbitrarily chosen twisted and untwisted quark masses. 
In order to make contact to numerical lattice data one has to properly match the untwisted quark 
mass according to the chosen definition for the critical quark mass. With the so-called parity 
conservation definition [21,22] one tunes m such that g>wt = which implies t = 0. In this case 
the expressions for the pion mass and decay constant reduce to what one expects from continuum 
ChPT With the pion mass definition one sets the bare untwisted mass to the value where the pion 
becomes massless (for vanishing pi). In this case the gap equation implies that t depends on a and 



therefore on pi. Hence the presence of t in ( [2.5D - ( |2.7[ ) causes deviations from the continuum ChPT 
expressions. 

Stalling from the gap equation one can derive some general characteristics of t(a). The ex- 
pected behavior is sketched in fig. [T[ For a — > one finds t — ► 1 , while for a — > °° it converges to 
the solution of the equation t — j8(l — 2? 2 )=0 which satisfies |?| < 1. The sign of ?_ is equal to 
the sign of j3. If j8 < there exists a value a where t = and &>wt = tc/2. Increasing the twisted 
mass beyond a results in angles &>wt > tt/2, and such values are indeed observed in the data [4]. 
We emphasize that without the O(pia) terms one finds the bound &>wt < 7r/2 for all a [9]. 
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Figure 1: Sketch of the solution t (a) with the pion definition for the critical mass. Depending on the sign 
of j3 the solution t converges to either a positive or negative value for large a. For j3 < there exists a finite 
a where t — and consequently 0)wt = ft/2, t (0) = 1 independently of the sign of j3. 

The 0(/j.a) contributions are also important for describing another numerical result. The parity 
conservation definition tunes mo such that t = 0. If we require t = to be a solution of the gap 
equation we find the condition 



m 



A4 



(2.8) 



This equation predicts a linear dependence between the bare masses niQ and which is also 
observed in the data [4]. 



3. Analyzing the data 

The data in Refs. [3-6] was generated with the Wilson plaquette action at /3 = 5.85 (a « 
0.123fm) and j3 = 6.0 (a as 0.093fm). The standard Wilson fermion action with a twisted mass 
term was used to compute the mesonic quantities. The values for the twisted quark mass cover 
the range of m n /m p 0.3 — 0.8. The untwisted bare quark mass was tuned in two different ways, 
according to the pion mass and the parity conservation definition for the critical quark mass. 

We performed combined fits to the data for the ratio R = a 2 m 2 n /apL, af % and cotft)wT- The fit 



functions (2.5) - (2.7) contain 10 free fit parameters. A total number of 34 data points are available 
for each a. We obtain good fit results with % 2 jd.o.f « 0.3 — 1.3 if we exclude the heaviest mass 
data. Fig. H| displays some fit results for j8 = 6.0. Three fits are shown for R and af % , depending on 
the fit range. Obviously, WChPT describes the bending for small quark masses quite well. 

The fit parameters are determined with moderate errors. The parameters in the gap equation, 
which essentially describe the bending, have statistical errors of 0(5%). A source of larger uncer- 
tainty is the choice for the fit range. Depending on data points included in the fit the central value 
for some parameters varies by about 20%. This suggests that higher order corrections in the chiral 
expansion are not negligible, even though the results (|2.5|) 



( |2.7[ ) are able to describe the data quite 
well. This is not unexpected, since the smallest /x max in fig. ^corresponds to m K /m p « 0.65, which 
is fairly heavy. 



4. Conclusions 

We interpret our preliminary results as follows. 
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Figure 2: Fit results for the j3 = 6.0 data. Upper left panel shows amo as a function of afiQ when tuned 
according to the parity conservation definition. The remaining panels show R — a 2 m^/ajx, af% and cotOwT 
for the pion definition. The three lines show the results for different fit ranges when data with jio < /i max is 
included in the fit (only one line is shown for cotOJwT since all three lines would lie on top of each other). 

WChPT is indeed able to describe the bending phenomenon observed in quenched tmLQCD. 
The origin of the bending is - in the language of the chiral effective theory - the contribution of 
the non-trivial ground state, encoded in the function t(/i). The terms of 0(/ia) cannot be ignored 
at the lattice spacings we considered. Once these terms are included WChPT explains a variety of 
features observed in the data: The angle G)wt can assume values larger than n/2, the untwisted 
mass depends linearly on the twisted mass when it is tuned according to the parity conservation 
definition, and the pion mass and decay constant show the strong bending for small quark masses. 

Fits to the data provide estimates for particular combinations of low-energy parameters with 
good statistical precision but sizable systematic errors due to the choice of fit range. The analysis 
suggests that higher order corrections might be necessary in order to obtain reliable estimates for 
these parameters with small errors. Our final results will be given elsewhere [23]. 
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